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Abstract 

A simple axion model is proposed in the scenario of large extra dimensions 
where gravity scale is as low as 1 TeV. To obtain an intermediate-scale decay 
constant of the axion, the axion is assumed to live in a sub-spacetime (brane) 
of the whole bulk. In this model there appear Kaluza-Klein modes of the axion 
which have stronger interaction than those of the graviton. The axion brane 
plays a role of absorber of the graviton Kaluza-Klein modes. Because of these 
reasons phenomenology and cosmology of the axion become highly non-trivial. 
We discuss various cosmological constraints as well as astrophysical ones and 
show that the model is viable for certain choices of the dimensionality of the 
axion brane. The structure of the model proposed here provides a viable 
realization of the fat brane idea to relax otherwise very severe cosmological 
constraints. 
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I. INTRODUCTION 



It has been suggested that the fundamental scale of nature can be as low as TeV, whereas 
the largeness of the effective Planck scale or the weakness of the gravity in a long distance 
can be explained by introducing large extra dimensions When there exist n of such 

extra dimensions, the relation between the gravitational constant SttG^v = l/^pi(4) in 4 
dimensions and the fundamental scale Mpi(4+„) in (4 + n) dimensions is given 

Km ~ Km;+2_,„), (1) 

where Vn is the volume of the extra dimensional space. For Mp;(4+„) ~ 1 TeV, the size of 
the extra dimensions r„ is computed as 

r„ ^ -1/1/" ^ 10^-^MTe|"^MeV"^ ~ 10^"^^MTe|"^cm, (2) 

where Mtev = Mp;(4+„)/TeV. The case n = 1 is excluded because the gravitational law 
would change at macroscopic level, but the cases n>2 are allowed by gravity experiments. 

There are also astrophysical and cosmological bounds because there are a tower of 
graviton Kaluza-Klein (KK) modes which contribute to the supernova cooling 0, the total 
mass of the universe and the late-time photon production |0 etc. We may avoid some of 
these problems by assuming an extremely low reheating temperature of the early universe. 

To have an inflation model with such a low temperature is already a tough subject, espe- 
cially with the cosmological moduli problem 0. Furthermore, we need baryon asymmetry 
after the inflation while the temperature is quite low and we knew that the process which 
generates baryon asymmetry should not induce proton decay. It is better to have higher 
reheating temperature for many reasons, but it is very hard to avoid known restrictions. 

There is another problem in TeV scale gravity. The simplest model of this type would not 
provide intermediate scales which are necessary to explain phenomenological issues like, the 
apparent gauge coupling unification , small neutrino masses and the strong CP problem. It 
was pointed out 0] that particles dwelling in the extra dimensions, other than graviton, can 
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have similar effective interaction terms in four-dimensional physics. This type of interaction 
is suppressed as gravity and thus give extremely weak interaction between bulk matter and 
normal matters. For the neutrino, there have been many discussions and researches about 
how to extend the minimal setting to obtain the small neutrino masses and mixing p|-p!T|, 
but the axion as a solution of the strong CP solution has not been studied thoroughly in 
this context. In this paper we will propose an axion model in TeV scale gravity for various 
numbers of the extra dimensions, with special emphasis on cosmological constraints: the 
addition of the axion in the large-extra-dimension model substantially alters the cosmology. 

In the following sections, we will use the convention Mpi = Mpn^^) = reduced Planck mass 
and = Mp;(4+„) = fundamental scale. 



II. PQ SCALE IN EXTRA-DIMENSION 

If an axion is a boundary field confined on 4 dimensions, the Peccei-Quinn (PQ) scale 
fpQ is bounded by ~ 1 TeV. To obtain a higher PQ scale, the axion field has to be 
inevitably a bulk field. If it lives in the whole (4+n) dimensional bulk, the PQ scale will be 

fpQ ~ Mpl. 

However, the damped coherent oscillation of the axion vacuum at the early universe 
with fpQ ~ Mpl would overdose the universe. A conventional argument gives an upper 
bound fpQ < 10^^ GeV. Even when an entropy production takes place after the QCD phase 
transition {e.g. the reheating temperature is smaller than ~ 1 GeV), fpq cannot be much 
larger than 10^^ GeV. 



In Ref. [0, it was shown that the late thermal inflation can raise bound of the PQ 
scale up to 10^^ GeV. This argument can be applied to our case, provided that the coherent 
oscillation of the infiaton is followed by the reheating process. Then the bound on fpq is 

/.«<10"Oev(A)(^)(-)- 

where h is from Hubble constant in units of 100km sec~^Mpc~^, 6 is the initial value of PQ 
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vacuum angle and Tr is the reheating temperature after inflation. Thus the axion model 
suffers from the over-closure problem if there is no other scale than and Mpi. 

In this paper, we propose that a natural way to realize an intermediate scale axion is to 
make it live in a (4 + m) dimensional sub-spacetime ((3 + m) brane) (m < n) of the whole 
(4 + n) dimensional bulk. The idea of using sub-spacetime to realize an intermediate scale 



already appeared in Ref. |I0| in the context of neutrinos. 

Let X be a complex scalar field which contains PQ axion in 4+m dimension; a. If the 
axion field lives only on 4 + m dimensional sub-spacetime where m < n and the volume of 
extra-dimension is Kn, 

C^ = j dx^+^'d'"x*dMX + / (4) 

where means an extra-dimension coordinate. Assuming that the vacuum expectation 
value of the x fi^ld does not depend on the extra-dimension coordinates, we obtain 

^ 103(1+5™/-) M^;™/"GeV. (5) 



Here we have defined the 4-D axion field as a = ^/Vrnfl-ix^ = 0) and assumed that the size 
r„ is common for all extra dimensions. 

A lower bound of fpQ comes from astrophysical observations, e.g. red giant and super- 
nova cooling by axion emission. It is known that fpq should be larger than lO^GeV. In 
extra dimension physics, KK modes also contribute to supernova cooling if their masses are 
smaller than the core temperature (~ 30MeV). We discuss it in the next section. 

To have 10^ GeV < fpq < 10^^ GeV, we need 2/5 < m/n < 4/5. Possible sets of {m,n) 
with fpQ where M^, = 1 TeV and 10 TeV can be found in Table |, ||. 



III. LABORATORY AND ASTROPHYSICAL CONSTRAINTS 

Experiments on detecting axion from the nuclear reactor and the sun give bounds on 
PQ scale < 10^ GeV (Lab), 10^ GeV < fpq < 10^ GeV from the sun. However this limit 
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strongly depends on the photon axion couphng and method of detecting axion. Furthermore 
the laboratory bound is much weaker than the astrophysical bound. 

The strongest bound from astrophysics is the Supernova cooling. In SN1987A observa- 
tion, it was calculated that 

-J- < GeV-l (6) 

JPQ 

Since the axion KK modes interact exactly the same way as the conventional axion, the 
effective interaction of the KK modes at the core temperature T ~ 30 MeV is 

1 T^m 

1 1 18 r^„^.-2 



f2 X (^--r - < 10-^^ GeV- (7) 

JPQ 

This gives a bound on the fundamental scale 

> (10^^ X 0.03™)^ GeV. (8) 

For m = 1, > 300 TeV [| and m = 2, > 5 TeV. ~ 1 TeV is allowed only if 
m > 2. This suggests that the number of extra-dimension should be at least 3 to include 
PQ mechanism as a solution of strong CP problem in TeV scale gravity model. 

In near future, high energy accelerator experiments may probe the axion KK mode 



emission. The graviton KK mode signals in collider were discussed in detail recently. |T4 
The scattering cross section of the graviton KK mode emission from high energy scattering 
with center of mass frame energy y/s is 



a oc 



1 / /~\ "+2 1 



(v^-'^)'^ -[tt] (9) 



'pi 

while the KK axion production is 



My-'" \MJ s 



1 / ^\ '"^^ 1 

, . ^,yi.„r . _. (10) 



^It was argued that there arise large quantum corrections in the case of one extra dimension in 



general, which would destabilize the gauge hierarchy |13|. This argument would exclude the m = 1 
case. 



Since the energy dependence of the axion KK mode cross-section is different from the gravi- 
ton KK mode cross-section, it might be possible to detect this difference at TeV scale collider 
experiments. 



IV. THERMAL PRODUCTION OF AXION KK MODE 



Since the axion dwells in the extra-dimensional brane, their masses are proportional to 
and they have stronger couplings than the graviton KK modes to the normal matters 
in general cases. 

If there is no hidden particle which couples to the axion or the mass of the KK mode 
axion is lower than the sum of three pion masses ~ 500 MeV, the main decay channel of 
light KK axion is to two photons. The life-time of each KK mode for a given (n, m) can be 
found in Table |, 0. 



coupling which is usually within 0.1 to 1. This decay can be cosmologically dangerous. For 
instance, for fpQ = 10^^ GeV and rriA = 1 MeV, life time of KK mode is ta ~ 10^'' sec, 
which is about the age of the universe. 

The graviton KK modes have similar cosmological problems because they can overdose 
our universe or decay into the photons at a late stage of cosmological evolution. Originally it 
was suggested that a "fat brane" can solve the cosmological problems by absorbing most 
of the decay products of the KK modes. However massless particles in the higher dimensional 
brane are not massless in our four-dimensional universe since they have momenta in the extra 
dimensions which appear masses in our universe 0. Another way to avoid this difficulty 
is assuming a large number of four-dimensional lattices in the bulk. But we need at least 
10^ empty universes. Or we should assume the existence of a four-dimensional hidden sector 
which has 10® times more degrees of freedom than those of the standard model while they 
should not be produced significantly by the reheating process after the inflation. 




(11) 



where is the mass of the axion KK mode and is the model dependent axion-photon 
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If we add an axion as a "brane particle" in this model, the graviton "bulk particle" will 
decay to the "brane" axion more efficiently, since its decay width will be enhanced by factor 
{Morn}"^- The graviton KK mode with mass mo will have decay width to the axion 

r{gKK ^ 2a) ^ O{10-')-§ x (mcr^)-. (12) 

After some period, instead of the massive graviton KK fields, we will have the axion 
KK fields with about same masses. Since the massive axion KK mode can decay into 
the photon pairs and the life time of the graviton KK mode becomes much shorter, the 
primordial graviton KK mode will not over-close the universe. Instead, it will contribute to 
the cosmological background radiation. This can be a severe constraint to the axion model. 
Also axions can be produced thermally during the reheating process. For these reasons, we 
have to check whether the axion model can survive the cosmological constraints. 

To estimate the constraints for various cases, we calculate the amount of the thermal 
axion produced at the reheating temperature and the axion KK modes from the graviton 
KK modes decay. In Appendix A, we derived the Boltzmann Eq. for the yield Y = p/s for 
four different sources of the axion KK modes: 

I. decay of graviton KK modes from the inverse decay, {2j, e'^e~, Dv — > Qkk)', the energy of 
the KK modes are concentrated on uia ~ Tr, 

II. decay of graviton KK modes produced from the scattering, (e7 — > egxK, e^e~ — > ^gKK)] 
this contribution is significant only if < T/j, 

III. the axion KK mode from the pion scattering (tttt — > naKx)', this process dominates if 
Tr > 10 MeV 

IV. the axion KK mode from the two photon inverse decay (27 — > ukk) ', this gives significant 
contribution when ~ Tr. 

For each case, 
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where 



^.=cj,i^r'7!^r. (20) 



3/2 

For the details of these calculations and the definition of function /(T), see appendix A (we 
present the numerical plot of I{T) in Fig. 5). 



V. COSMOLOGICAL CONSTRAINTS 

In this section we would like to discuss various cosmological constraints on the model for 
a given (m, n) set. Notice that in Table I and II, n > 5 in both M* = 1 and 10 TeV case 
and n = 4 in 1 TeV are cosmologically safe if Tr is low enough (~ 1 MeV). If the minimal 
KK mode mass is greater than 1 MeV, the KK modes are not generated in the thermal 
bath of such a low reheating temperature. On the other hand m < 2 in M* = 1 TeV and 
m = 1 in = 10 TeV is forbidden by the astrophysical bound. The cases n = 4, m = 3 at 

= 1 TeV and n = 3, m = 2 at = 10 TeV are not trivially allowed or ruled out by the 
cosmological constraints. Details on these cases are discussed in Appendix B. 



A. Big bang nucleosynthesis 



At the temperature of the universe around 1 MeV, it is required that there should not 
be additional particles which contribute to the energy density significantly. Otherwise ^He 
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would be produced more than what is observed now because the universe should expand 
faster than the standard scenario. 

We apply a rather loose bound that the energy density contribution by the KK mode 
should be smaller than one neutrino energy density at T = 1 MeV. At high Tr(> 10 MeV), 
axion KK mode production dominated by process III. In case III, one may practically have 
maximal mass 

mi = max{m7r, Tr} (21) 

for the KK mode which is produced in the thermal bath. 
Thus the bound from BBN is 



PA 

s 



~ D / dniA X (mArn)"^YA 

BBN Jmo 

m+1 



~ miF3(mirJ- ~ 2 • IQ-^^^^^A^ < 0.1 MeV, (22) 

where D is normalization constant which is equal to m in torus compactification with uni- 
versal distance r„. In our calculations we restrict ourselves to this case. 
If Tr ~ mjr, approximately 

20 

> 10^ X m^. (23) 

For m = 1, this reads > 600 TeV, and for m = 2, > 15 TeV. But if the reheating 
temperature is as low as 10 MeV, this bound is not important since As{Tr) is suppressed 
exponentially. 

B. Over-closure of Universe 

The total energy of the axion KK modes at present must not exceed the critical density: 

PA< Pc = ^- lO'^Soh^ MeV, (24) 
where sq — 3000cm~^ is the entropy of the present universe. 



For the case that the KK modes decay into some relativistic particles, we can divide the 
bound in two parts; decay before the present time and do not decay till now: 

^ ~ D p drriAimArnrYA + D H dmA{mArnT^f^ < 3 ■ 10~'h^ MeV, (25) 

So J mo Jm2 T[mA) 

where the axion KK mode with mass 1712 decays at the present time, and mi is defined in 
Eq. (|1D. 



C. Light element destruction 

Energetic photons from heavy axion KK modes which decay after 10^ sec can destroy 
the light elements made during the nucleosynthesis. Therefore there are several bounds on 



the density of the KK modes which weigh than 10 MeV. They are |15 



^ < lO-i^GeV (26) 

for > 10^ sec, 

— < 10"^ - 10-^°GeV (27) 
s 

for 10^ < ta < 10^ sec. Usually these bounds are less important than cosmological microwave 
background bound given below. 



D. Cosmological Microwave Background Radiation 

If the massive KK modes decay after ta > 10^ sec but before the recombination era, the 
produced photons may give a distortion of the cosmological microwave background radiation. 
The COBE observation gives a bound 

^ < 2.5 X 10-^Td (28) 



where Tq is temperature at KK mode decay. 
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E. Diffuse photon background 



Observations of diffuse photon backgrounds at the present universe give upper bounds 
on additional contributions to photon spectrum. For example, for the energy range 800 keV 
<E <30 MeV |T7[ 

^.Bx^(W"-78(^)"". (29) 



Constraints on other ranges of the photon energy can be found, e.g. in Ref. fTE 
Theoretical prediction is 

dJ^ uac , . 

3?T = ^ Br (30) 
dn An ^ ^ 

for the life-time of KK mode shorter than the age of the universe, and 



dj^ ^ nAcTa^^-.2^, ( 2E 

— Br X — 

dVt An Br Ho Xttlac^ 



Srx^i^h^ (,n,r.r (31) 



for its life time longer than the age of the universe |]T9|. Here we have introduced Br as a 
branching ratio of axion decay into two photons. 

In our brane picture, there is a priori no reason that the axion brane contains only 
our four-dimensional wall. Rather it will be natural to imagine that there is a parallel 
universe(s) or another four-dimensional wall in the brane. Or one can just imagine that 
there are some unknown particles on our wall itself. Then one can consider the situation 
that thermally produced KK modes of graviton will mainly decay into the axion in the brane 
and axions both from the graviton decay and the thermal production will decay into the 
parallel wall if this wall has some kind of QCD and/or f/(l) type of interactions. (Or it can 
decays into some hidden QCD/f/(l) fields in our universe.) Since the axion decay width 
is highly suppressed by {^^^ with a being the fine structure constant, it is easy to get a 
low branching ratio to decay into the photon, in other words, most of the axions decay into 
the other wall(invisible section) if the coupling constant, the color factor of the other gauge 
interaction, and/or the number of the fermions with PQ charges in the decay loop diagram 
in the other wall(invisible section) are large enough. 
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F. Results 



Here we summarize the results we obtained. The reader should be referred to Appendix 
B for more detail. (We approximate h = 0.7) 
1. n — A, m — 3 and = 1 TeV case: 

1. BBN bound 

Tr < 90 A3 MeV ~ 80 MeV. (32) 

2. Over-closure bound, 
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Tr < 12 (^-^^j MeV ~ 30 MeV (for Br^l). (33) 

3. CMBR bound, (for tua > 100 MeV and Tr > 10 MeV) 

Tr<2- 10-' ' MeV. (34) 

4. Diffused photon bound, 

for Tr > 10 MeV, 

-2 / BrY^'^^ ( \"°-^^ .-1/3 



for Tr < 10 MeV, 

Tr < 0.3Sr-°-^^ MeV, (36) 

where Br < ra-^2j/Ho, and 

r,i<5C-°-^«MeV (37) 

where Br > Va^2T/HQ. 

II. n = 3, m = 2 and = 10 TeV case: 

1. BBN bound 

Tr < 100 A~^^^ MeV ~ 90 MeV. (38) 
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2. Over-closure bound, 



rR<28(-^|^j MeV~40MeV (for 5r = 1). (39) 



3. CMBR bound, 



Tr<^- 10-^ ' MeV. (40) 



4. Diffused pfioton bound, 
for Tr > 10 MeV, 



T« < 3 . 10- f ^) f ^-^^ MeV, (41) 

" lC„^/ VlOMeV; 3 , V ; 



0.1 , 



for Tr < 10 MeV, 



where Br < ra_>2-y/i?0) and 



< O.lEr"'-^ MeV, (42) 



Tr < 3 MeV (43) 



where Br > Ta^2'y/Ho. 

In the both cases, n — A, m — 3 and M* = 1 TeV and n = 3, m — 2 and = 10 TeV, 
hfe times of axion KK modes up to ~ 1 GeV are quite long and so the constraint from 
the light element destruction will not give any further bound. 

We also performed computer calculations on both n — A, m — 3 and — 1 TeV and 
n — 3, m — 2 and = 10 TeV cases. Fig. 1 and 2 show that the DPBR is the stringent 
bound if the branching ratio is small. But in an extremely small branching ratio case, CMBR 
is dominant. This is because the life time of the axion KK modes becomes shorter than 10^^ 
sec, and so the produced photons will disturb the CMBR spectrum. One can see that BBN 
is independent of the branching ratio. Note that these behaviors are consistent with what 
we observed from the analytical computations given above. 

In Fig. 3 and 4, we present combined cosmological limits for all possible combinations of 
n, m in Table 1 and 11. We find that the reheating temperature is allowed to be significantly 
large if the branching ratio to the photon is small enough. 
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VI. CONCLUSIONS 



In this paper, we discussed the axion model in the extra dimensions whose PQ scale 
lies in an intermediate scale fpQ < 10^^ GeV. This intermediate scale can be obtained by 
introducing a 3 + m dimensional brane in the 4 + n dimension bulk. 

If we include the axion as a brane particle, it will change the phenomenology of the 
extra dimension physics, especially, cosmology. Since the graviton KK mode will decay into 
the axion KK mode, the over-closure problem is not as serious as the original model of 
Arkani-Hamed et.al. On the other hand, the argument from stars and supernova cooling 
will give a more strict bound on the axion production. Among other things, the most severe 
cosmological bound comes from photon emission through the decays of the KK modes of the 
axion. We found that the astrophysical argument restrict the number of the dimensionality 
of the sub-spacetime where the axion lives: m > 2 for = 1 TeV and m > 1 for = 10 
TeV. The latter cosmological argument requires quite a low reheating temperature after the 
inflation. 

To lift this bound, we can introduce the hidden matter/gauge fields to another four- 
dimensional wall (or even to our wall itself) which has much stronger coupling to axion 
and/or much more generations of particles, (or maybe much lower QCD phase transition 
scale, etc). This can significantly lower the branch ratio of the axion KK mode decay into 
photons. 

The whole picture can be used to improve the original fat brane idea. This higher 
dimensional object plays a role of an absorber of the KK graviton modes. If the fat brane 
couples to four dimensional wall(s) with interaction stronger than gravity, the produced 
particles in the fat brane may then decay into relativistic particles on the four dimensional 
wall(s). This mechanism can solve the problem of the overclosure of the universe by the KK 
modes. Note that the fat-brane particles are not necessarily the axions as we discussed, but 
can be any other weakly interacting particles living in a higher dimensional brane. Moreover 
if the particles produced by the fat-brane particle decays do not contain photons or any 
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other cosmologically dangerous particles, then we can avoid other cosmological problems 
such as the ones related to the cosmic photon backgrounds. It is worth mentioning that 
some mechanisms of generating small neutrino masses have a similar structure, which will 



be discussed elsewhere 20 



In this paper, we did not consider an alternative solution to the strong CP problem such 
as a spontaneous CP breaking model. This typically requires very heavy quarks. If we 
want to keep all particles carrying standard-model gauge charges in four dimensions, the 
maximum scale we can have is the fundamental scale i.e. around TeV. For this reason, this 
type of model is not favorable in the extra-dimension scheme. 
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APPENDIX A: ESTIMATION OF THE PRIMORDIAL AXION KK MODE 

DENSITY 

The production rate of the axion KK mode from initial particle i and j can be calculated 
from Boltzmann Eq., 

riA + 3H{t)nA = '^{(Tv)ijninj, (Al) 

where n^, H{t), a, rii represent axion KK mode number density, Hubble constant, scattering 
cross section and number density of initial particle i. 

Using the relation t = 0.5H{T)^^, we can convert the time parameter to inverse temper- 

(\ 1/2 rp2 
Y^j jf-^. If we assume that the KK modes are produced 
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by the particles in equilibrium, we can rewrite ( [Al| ) with the yield Y = ua/s. {s is entropy 
density s = ^g^sT^, where (7*5 ~ is approximately 10 for 1 MeV < T < 100 MeV.) 

dY X 



dx H{m 

where 



-TaYeq (A2) 



^A = nEQj2i'^^)ij- (^3) 



A similar equation can be derived for the inverse decay case. 

The yield of the KK mode at the equilibrium Yeq is about 0.28/(?^,s when initial particles 
are relativistic (or proportional to exp(— x) if they are non-relativistic). After integrating 
Eq. ([A^) and ( |A4| ) from the reheating temperature Tr to present temperature, we will get 
result with form 

which can be used most of calculations reliably. 

Let's estimate the sources of the axion KK mode production. The KK modes of axion 
can be produced from either the thermal graviton KK mode decay or initial thermal bath. 
We will classify four relevant cases. 

Class I: The KK mode of graviton which has mass around the reheating temperature 
me — Tji generated dominantly through the inverse decay 77 — > gxK, ^'^ Qkk and 
e+e" gxK 

= (A6) 

__ ml 

for initial spin averaged. For the tensor mode of graviton KK mode we should multiply 5 to 
r. This will generate 
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where 



Here we approximated that ~ after the graviton decay. This mode is most abundant 
at rriA = Tr and decrease fast if ^ Tr. 

Class II: The KK modes which has much less mass than mc < Tr will be produced 
dominantly by the scattering processes e^7 Qkk^^^ e^e~ — > QkrI- If we choose a limit 
that KK mode mass is less than Tr but greater than the electron mass, we can calculate 



the interaction rate from the amplitude presented in [|4 



r.+.--»„, --^J- (('■3)J'r (All) 



The yield from scattering is 



¥2-2- 10-23 ( \ A2, (A12) 

VlOOMeV/ ' ^ ' 



where 



for each KK mode with mass itia- This bound is valid only if A2 is positive, i.e. the reheating 
temperature is significantly higher than the KK mode mass. 

Class III: Thermal axion produced mainly by pion-pion scattering 7r^7r° — > vr^a, vt+tt^ — >■ 
vrOa, for Tr > 10 MeV. 

r,„„y,, . (^) X i(T) (AM) 

where C^^ is (1 — z)/3{l + with z = mu/md. In the limit of = 0, we can used the 



temperature dependent function I{T) in Ref. ^1 
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/(T) = / dx,dx,^fiy,)fiy,) j" ^^i^^I^l^^jKl. (A15) 

J yiy2 J-i s^T'^ 

where f{y) = l/(e^-l), = bH/T, y^ = E^/T {i = 1,2) and s = 2{ml + T\y^y2-XiX2Uj)). 
We justify ourselves using tua — limit by presenting the plot produced by computer (Fig. 
5) which shows that the mass dependence of I{T) is indeed small. 

I{T) is around 10^ for T > 50 MeV, and it suppressed exponentially at T < 10 MeV. 
For T <C m^, we can approximate function /(T), 



8 \ T J ^ \ T 
We can estimate thermal axion KK mode yield from pion scattering, 



where 



Class IV: Thermal axion can be generated through the two photon inverse decay, 

r , r^^f^VnA ^27- lo-'C ^ (Am 



which leads 



PpQ \ fpQ ) VlOOMeV/ 

where 



This shows that the axion produced from the inverse decay will dominate over the axion 
from the pion scattering III, if Tr < 10 MeV. Thermal axion can be produce with photon 
electron scattering process II, we will ignore this axion unless the reheating temperature is 
very low (compared with the process IV, it should be less than a few MeV). 



To estimate the energy density of the axion KK mode for given T, multiply the entropy 
density at the temperature s(T) and the mass of axion KK mode to the yield. We should 
count the number of KK modes for the allowed energy range, (£^r„)^ where E is typically 
Tr. But for case III, Tr — max{m^, Tr}. N is the dimension of the bulk where the produced 
particle exists. N — n and N — m ioi case {1, 11} and {111, IV} respectively. 

APPENDIX B: CALCULATION OF THE COSMOLOGICAL BOUNDS 

In this section we describe the details of these calculations for two non-trivial cases, 
n = 4, m = 3 at = 1 TeV and n = 3, m = 2 at — 10 TeV, though we have calculated 
with computer all relevant cosmological bounds with allowed sets of m and n and presented 
it in Fig. 3 and 4. 

This is the case with relatively low Tr < C(IOO) MeV, because of large fpg ~ 10^^ 
GeV. Since the life-time of the axion KK mode is quite long for ttia < 100 MeV and CMBR 
is stronger bound than the light element breaking bound, the bound from light element 
breaking is not relevant in these cases. 

1. Big Bang Nucleosynthesis Bound 

The axion KK mode with mass below 100 MeV cannot decay before 1 sec and will be 
restricted by BBN bound on neutrino species. Let's assume that Tr > lOMeV then 



Pa 
s 



BBN J mo 



rmi 

D / dm A (mArn)"'YA 

J mn 



^ m,Y^{m,rr,r ~ 2 • lQ-^^-^h_A^ < 0.1 MeV (Bl) 



fE2 

D / dm A X = # of KK modes between Ei and E2. 

J El 



where mi = max{m7r, Tr} and D is normalization constant defined as 

rE2 

I El 

We used the torus compactification with uniform distance approximation Vn = r^, Vm = T'n 
so that D — m. This gives a bound for n = 4, m = 3 at 1 TeV 
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Tr < 90 A^^'^ MeV ~ 80 MeV (B2) 
and for n = 3, m = 2 at 10 TeV 

Tr < 100 A^^'^ MeV - 90 MeV. (B3) 

2. Over-closure Bound 

We can divide this bound by KK mass 1712 ■ If > ^2-, its life-time is shorter than 
the age of universe, otherwise it will remains as a cold dark matter. For Tr < 10 MeV, the 
over-closure is not a problem in the region we are interested. Therefore Y3 is most dominant 
source of axion KK mode in this case. The total cold dark matter density will be 



^ ~ T) dniA {niArnTY^ < 3 ■ lO^^/i^ MeV (B4) 



where m2 — 10 MeV, mg = ^ 30 keV for = 4, m = 3 at 1 TeV and m2 — 1 MeV, 
mo ~ 5 keV for n = 3, m = 2 at 10 TeV. Using Eq.( [XT7D , Eq.(|B|) becomes 

2 ■ 10-3m^+^— ^^^3 < 3 • IQ-^h^ MeV. (B5) 

J TT * 

We can get the bound on Tr for n = 4, m = 3 at 1 TeV, 

/ \2/3 

Tr < 70 { — ] A3 MeV (B6) 



0.7 



and for n = 3, m = 2 at 10 TeV 



Tr < 350 i—j A,'/'MeY. (B7) 

We can consider another situation that the axion KK mode which can decay into some 
relativistic dark matter X dominantly. In this case, 

^ ~ D / dniA {mATnTY^TFT-^ < 3 " 10"'/^' MeV, (B8) 
So Jm2 T{mA) 

where Tq ~ 2 x 10~^^ GeV is current temperature of universe and 
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T(m^) ^ JVa^^^xMpi ~ 2 ■ 10-4a^5r-i/2— i-^ (B9) 



is the temperature when axion KK mode with mass decays, where Br is the branching 
ratio for 

which is the case that majority of axion KK modes decays into the invisible relativistic dark 



matter X. Then Eg . ( p8|) becomes 

2 • 10"^5r^/2ZfQ ^3 MeV < 3 • lO^^/i^ MeV. (Bll) 

For the reasonable range of Br, we can set mi = rriT,. For n = 4, m = 3 at 1 TeV, 

f h 

Tr < l2Br-^'^ { — j Cli^^A-^^'^ MeV ~ 30 MeV (for 5r = 1), (B12) 



and for n = 3, m = 2 at 10 TeV, 



/ \2/3 

Tr < 28Br'^/^ i — j MeV ~ 40 MeV (for Br = l). (B13) 

3. Cosmological Microwave Background Radiation 

CMBR bound is the most severe during the time period 10® < < 10^^ sec. The bound 

^ < 2.5 X 10-^T^ (B14) 



is actually weaker than other constraints if the branching ratio Br is large. In this case, the 
reheating temperature should be relatively small Tr < 10 MeV. But if Br is very small, it 
gives stronger bound than other cosmological constraints. is the same as Eg. (|B9|) . If we 
assume Tr > 10 MeV, and set Ap^ ~ Br x pa(Td), Eg.( pi4|) becomes 

2 . lO-^Er . mT'^^,A, < 5 ■ lO'^^.i^r-V^ ^^V (B15) 
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If Tr < mTT, the maximal value of KK mode mass is around m-j^. This leads the bound, for 
n = 4, m = 3 at 1 TeV, 

Tr<2- lO-^Sr-^/^Ag ^/^C„^/=^ MeV (B16) 

for n = 3, m = 2 at 10 TeV, 

Tr < 4 • 10-^ Br-^^^As^^^Clf^^ MeV. (B17) 

To have Tr > 100 MeV, we need a very small branching ratio ~ 10~^. (This bound is not 
valid if Tr < 10 MeV.) 



4. Diffused Photon Background 

Since the hfe-time of axion KK mode is longer than 10^^ sec in allowed region, we 
obtain most strong bound on reheating temperature from diffused photon background. Let's 
consider three cases. 

Case 1. Tr > 10 MeV : In this case, the majority of KK mode will decay before present 
time. The observed bound when 800 keV < < 30 MeV is, 

1^<78(_±_) cm-^-W^ (B18) 
dfl VlkeV; ^ ' 

corresponds for theoretical prediction 

- = ^ i?r^(m^r„) . (B19) 

Here SqC ~ 9 • 10^^cm~^sec~^. This will give the inequality, 

Br xYs{mAr„r <6- 10-'' (^-^^ (B20) 

where E ~ 771^/(2(1 + z)) 

1 + z ~ 4 X 10^^ Uloh^y'^' fl£.y^' , (B21) 

The time of axion KK mode decay into two photons is 
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m^Br- r,-;^^2, ^ 3.7 ■ 10' Br^. (B22) 



''A 

Thus the present energy of diffused photon (for fio — 1) is 



Inserting this result into ( [B20| ) and setting h = 0.7 leads 



and it gives a bound for n = 4, m = 3 at 1 TeV, 



Tn<2. 10- ( ^) (-^\ A- MeV (B25) 



Ca^l VlOMeV/ 



and for n = 3, m = 2 at 10 TeV, 



UOMeV. 



Tr < 3 • 10-' — TTT^TTT ^3 MeV. (B26) 



Case 2. Tr < 10 MeV, Br < Ta^^^^jH^- 

In this case y ~ 14 and axion KK mode will decay before present time, i^o — 2 • 10~^'/i 
GeV is Hubble constant of present universe. Then the relation with energy and life time 
become same as Eq.( P24|) , 



<10io(^l C]^^Br-'-'(-^^X' (B27) 



^GeVy VlOMeV / 

and it gives a bound for ?T, = 4,m = 3at 1 TeV 

and in tjia ~ Tr limit, 

Tr < 0.3 5r~°-^3 MeV. (B29) 
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For n = 3, m = 2 at 10 TeV 



Br \ ^'^ / tua \ 



and in ttia ~ ^i? limit, 

T/j < 0.1 Br-^-^ MeV. (B31) 
Case 3. Tr < 10 MeV, Br > Ta^^^2j/Ho: 

In this case, axion KK modes will not decay before the present time. Then the theoretical 



prediction becomes (see references, for instance. Sec. 5.5 in ||19|| ) 



a\l An BrHo \mA/ 



(B32) 



and this will lead 



Approximately, for n = 4, m = 3 at 1 TeV, 

T« ( <2-10'(—] C-M,' ( -^V'" MeV. (B34) 
^VMeVy \0.7j ^ VMeVy ^ ^ 

We can approximate niA ^ Tr E, then ^4 ~ and 



h 



0.12 

-0.48 



Tr < 5.5 ^— j MeV. (B35) 

The case n = 3, m = 2 at 10 TeV, the majority of thermal axion KK mode with Tr > 1 
MeV will decay before present time. The condition Br ~ T aji^^2'y/ Hq and Eg. ( [B3 1|) will 
determine the reheating temperature 

Tr<2^?> MeV. (B36) 

We did full calculations for all relevant sets of m and n with computer and get consistent 
results. For instance, you can see three different regions of DPBR bound in Fig. 1 and 2. 
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TABLES 

TABLE I. PQ scales and lifetimes of axion and graviton KK modes for M* = 1 TeV, where 

— rriA 
100 — looMeV 



(m, n) 


r-i[MeV] 


/pQ[GeV] 


rp[sec] 


Ta [sec] 


(1,2) 


4 X 10-^° 


5 X 10^° 


2 X 10^M{^^ 


6 X lO^M^-^g 


(2,3) 


6 X 10-5 


2 X 10^3 


1 X lO^M^^^ 


8 X lO^^M^o^ 


(2,4) 


2 X 10-2 


5 X 10^° 


2 X lO^^Mf^,^ 


6 X 10'^ M{qI 


(3,4) 




3 X 10^^ 


3 X 10^M{^g 


3 X 10^^M{qI 


(2,5) 


0.7 


1 X 10^ 


2 X lO^^M^o^ 


5 X lO^M^^g 


(3,5) 




2 X 10^2 


1 X lO^^j^^fj^^ 


7 X 10i°M{^g 


(4,5) 




2 X 10^5 


9 X lO^Mf;,^ 


1 X 10i^M{o^ 


(3,6) 


7 


5 X 10^° 


2 X 10^5 Mio^ 


6 X 10^M{^^ 


(4,6) 




2 X 10^^ 


1 X lO^^Mi^^ 


8 X lO^^M^f^ 



TABLE IL PQ scales and lifetimes of axion and graviton KK modes for = 10 TeV. 



(m. ii) 






r,Jsoc] 


T„ [sec] 


(1,2) 


4 X 10-^ 


2 X 10^^ 


2 X lO^M^Q^ 


6 X W^M^qI 


(2,3) 


3 X 10-^ 


4 X 10^3 


3 X lO^Mi^^ 


4 X 10^3^100 


(2,4) 


0.6 


2 X 10^1 


2 X lO^'^M^f^l 


6 X lO^M^^g 


(3,4) 




6 X 10^^ 


1 X lO^^j^^f;,^ 


9 X IO^^M^qI 


(2,5) 


20 


6 X 10^ 


1 X lO^^M^o^ 


8 X lO^M^o^ 


(3,5) 




4 X 10^2 


2 X lO^^Mf^^ 


5 X lO^^Mf^g 


(4,5) 




3 X 10^5 


4 X lO^^Mf^^ 


3 X 10^^ M{qI 


(3,6) 


160 


2 X 10^1 


2 X lO^^M^o^ 


6 X 10^ Mfog 


(4,6) 




4 X 10^^ 


3 X lO^^Mi^^ 


4 X 10i=^M{^g 
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FIG. 1. The bound of Tr and Br for M* = 1 TeV and n = 4, m = 3, The upper and right side 
of each hne is excluded region. 
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FIG. 2. The bound of Tr and Br for = 10 TeV and n = 3, m = 2. The excluded region is 
same as Fig. 1. 
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FIG. 3. The bound of Tr and Br for M* = 1 TeV. The upper and right side of each hne is 
excluded region. 
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FIG. 4. The bound of Tr and Br for M* = 10 TeV. The excluded region is same as Fig. 3. 
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FIG. 5. The function I{T) for various ruA in the KK mode yield III. 
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